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The  detailed  analysis  of  measured  interferograms  generally  requires  phase  correction.  Phase-shift 
correction  methods  are  commonly  used  and  well  documented  for  conventional  Fourier-transform  spec¬ 
troscopy.  However,  measured  interferograms  can  show  additional  phase  errors,  depending  on  the 
optical  path  difference  and  signal  frequency,  which  we  call  phase  distortion.  In  spatial  heterodyne 
spectroscopy  they  can  be  caused,  for  instance,  by  optical  defects  or  image  distortions,  making  them  a 
characteristic  of  the  individual  spectrometer.  They  can  generally  be  corrected  without  significant  loss 
of  the  signal-to-noise  ratio.  We  present  a  technique  to  measure  phase  distortion  by  using  a  measured 
example  interferogram.  We  also  describe  a  technique  to  correct  for  phase  distortion  and  test  its  per¬ 
formance  by  using  a  simulation  with  a  near-UV  solar  spectrum.  We  find  that  for  our  measured  example 
interferogram  the  phase  distortion  is  small  and  nearly  frequency  independent.  Furthermore,  we  show 
that  the  presented  phase-correction  technique  is  especially  effective  for  apodized  interferograms.  ©  2004 
Optical  Society  of  America 
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1.  Introduction 

Spatial  heterodyne  spectroscopy  (SHS)  is  a  relatively 
novel  concept  that,  when  compared  with  other  spec¬ 
troscopic  techniques  such  as  Fabry-Perot  or  Michel- 
son  interferometers,  can  offer  many  advantages  for 
high-spectral-resolution  diffuse-source  spectroscopy. 
It  was  conceived  in  the  late  1980’s  and  was  made 
possible  primarily  by  the  availability  of  detector  ar¬ 
rays,  e.g.,  the  charge-coupled  device  (CCD),  in  com¬ 
bination  with  high  computing  speed  to  process  the 
recorded  interferogram  data.1 

The  basic  SHS  configuration  is  similar  to  a  Mich- 
elson  interferometer  with  the  return  mirrors  re¬ 
placed  by  fixed  diffraction  gratings.  It  contains  no 
moving  parts  and  can  be  field  widened  with  fixed 
prisms  in  the  interferometer  arms  to  increase  the 
etendue  of  the  instrument  by  as  much  as  2  orders  of 
magnitude  over  conventional  Fabry-Perot  or  scan- 
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ning  Michelson  interferometers.  Detailed  descrip¬ 
tions  of  SHS  can  be  found  in  publications  by 
Harlander  et  aZ.1-3 

Figure  1  shows  the  basic,  non-field-widened  SHS 
configuration.  A  major  difference  between  conven¬ 
tional  Fourier-transform  spectrometers  and  SHS  is 
that  in  SHS  all  interferogram  points  are  recorded 
simultaneously  by  an  array  of  detectors,  whereas  a 
Fourier-transform  spectrometer  records  the  inter- 
ferogram  samples  with  a  single  detector  in  a  time 
series  during  which  the  scan  mirror  is  moving.  In 
a  Fourier-transform  spectrometer,  light  from  the 
entire  aperture  of  the  interferometer  is  focused  on  a 
single  detector  element.  This  drives  the  need  for 
flat,  homogenous  interferometer  optics  and  a  well- 
controlled  scanning  mechanism  because  variations 
in  optical  path  difference  across  the  interferometer 
aperture  of  more  than  a  fraction  of  a  wavelength 
greatly  reduce  the  detected  fringe  contrast  and  the 
signal-to-noise  ratio  of  the  spectrum.  In  the  SHS 
configuration,  the  fringe  localization  plane  near  the 
gratings  is  focused  onto  the  imaging  detector,  as 
illustrated  for  one  detector  element  in  Fig.  1.  Con¬ 
sequently,  the  signal  for  each  interferogram  ele¬ 
ment  uses  only  a  small  area  of  the  gratings,  beam 
splitter,  and  the  exit  optics  before  it  reaches  the 
detector.  Variations  in  optical  flatness  across  the 
aperture  result  primarily  in  a  phase  distortion  in 
the  interferogram  and  produce  fringes  that  are  not 
straight  and  equally  spaced;  however,  the  fringes 
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Fig.  1.  Schematic  diagram  of  the  basic  non-field-widened  SHS 
configuration.  The  dashed  lines  illustrate  an  incoming  wave  front 
and  the  corresponding  exiting  wave  fronts,  which  are  crossed  with 
an  angle  of  2-y.  The  ray  bundle  for  one  interferogram  element  is 
outlined,  showing  that  only  a  small  section  of  the  interferometer 
and  optics  are  used  for  any  individual  interferogram  element. 
FTS,  Fourier-transform  spectrometer. 

are  recorded  with  high  contrast  as  long  as  the  in¬ 
terferometer  elements  are  good  over  the  small  area 
sampled  by  each  detector  pixel.  These  phase  dis¬ 
tortions  are  analogous  to  nonuniform  path  differ¬ 
ence  sampling  in  Fourier-transform  spectroscopy 
(FTS).  Another  potential  cause  for  phase  distor¬ 
tion  in  a  SHS  interferogram  is  the  image  distortion 
of  the  exit  optics,  which  reimages  the  fringe  local¬ 
ization  plane  on  the  detector  array.  These  phase 
distortions  have  no  obvious  equivalent  in  FTS. 
Generally,  significant  phase  distortion  can  and 
should  be  avoided  in  any  SHS  instrument  design. 
If,  however,  an  instrument  shows  significant  phase 
distortion,4  it  is  a  characteristic  of  the  spectrometer 
and  can  generally  be  corrected  without  significant 
loss  of  the  signal-to-noise  ratio  in  the  measured 
data  during  the  postprocessing  of  the  interfero¬ 
gram.  The  SHS  concept  therefore  poses  relaxed 
flatness,  index  homogeneity,  and  alignment  re¬ 
quirements  compared  with  a  Fourier-transform 
spectrometer.  Moreover,  optical  defects,  for  exam¬ 
ple,  a  scratch  on  the  grating,  have  only  a  localized 
effect  on  the  interferogram,  which  allows  for  their 
isolation  and  proper  treatment  in  the  processing  of 
the  interferogram. 

In  the  following  section  we  revisit  the  basics  of  the 
SHS  interferogram  as  a  basis  for  the  subsequent  dis¬ 
cussion  on  how  to  determine  the  phase  distortion  for 
a  given  spectrometer.  As  an  example,  an  interfero¬ 
gram  measured  with  the  Spatial  Heterodyne  Imager 
for  Mesospheric  Radicals-Space  Shuttle  Middeck 
(SHIMMER-MIDDECK)  instrument3  is  used.  Fi¬ 
nally,  correction  techniques  for  different  types  of 
phase  distortion  are  discussed  by  use  of  model  calcu¬ 
lations  of  a  near-UV  solar  spectrum. 


2.  Spatial  Heterodyne  Spectroscopy  Interferogram 

Figure  1  illustrates  the  basic  process  by  which  a  SHS 
interferogram  is  created.  Incoming  wave  fronts  are 
split  into  the  two  interferometer  arms,  and  the  re¬ 
turning  diffracted  wave  fronts  recombine  at  the  beam 
splitter.  The  recombining  wave  fronts  are  parallel 
and  coincident  for  only  one  wave  number,  the  Littrow 
wave  number  (u0),  which  produces  a  constant  signal 
across  the  detector.  For  wave  numbers  different 
from  the  Littrow  wave  number,  the  exiting  wave 
fronts  are  crossed  as  shown  in  Fig.  1,  resulting  in  a 
fringe  pattern  that  is  recorded  by  the  detector  array. 
For  wave  numbers  close  to  the  Littrow  wave  number, 
the  spatial  frequency  of  the  fringes  is  directly  propor¬ 
tional  to  the  wave-number  difference  (8u  =  <t0  —  ct),  so 
that  the  modulated  portion  of  the  detected  signal  is 
the  Fourier  transform  of  the  incident  spectral  distri¬ 
bution  heterodyned  around  the  Littrow  wave  num¬ 
ber.1 

For  the  following  discussion  of  phase  distortion,  we 
introduce  the  fundamental  form  of  a  SHS  interfero¬ 
gram,  in  which  x  is  the  location  of  an  individual  de¬ 
tector  element  in  the  dispersion  plane  [x  =  0  is  the 
point  of  stationary  phase,  i.e.,  the  location  of  zero 
path  difference  (ZPD)],  k  is  proportional  to  the  het¬ 
erodyned  wave  number,  8cr,  and  <L(k,  x)  is  an  additive 
phase-distortion  term: 


I(x)  =  J  A(k,  x)(3(K)(exp{j[i<x  +  d^K,  x)]} 

+  exp{— j[kx  +  <I>(k,  x)]})dK.  (1) 

The  goal  is  to  retrieve  the  spectral  intensity  (3(k)  from 
the  interferogram,  which  is  the  sum  of  cosine  fringes 
with  the  amplitude  (3(k)  multiplied  by  a  normalized 
envelope  function  A(k,  x).  If  the  envelope  function  is 
unity  and  the  phase-distortion  term,  <I>(k,  x),  is  zero, 
a  simple  Fourier  transform  yields  (3(k)  for  k  >  0  and 
f3( — k)  for  k  0.  The  effect  of  the  envelope  function 
is  that  the  retrieved  spectrum  is  convoluted  with  the 
Fourier  transform  of  the  envelope  function.  This  ef¬ 
fect  is  familiar  from  FTS  and  contributes  to  the  finite 
resolution  of  measured  spectra  and  the  instrumental 
line-shape  function  in  general.  In  SHS  the  envelope 
function  depends  on  wave-number-independent  con¬ 
tributions  such  as  the  finite  detector  size  and  wave¬ 
number-dependent  contributions  such  as  the 
modulation  transfer  function  of  the  imaging  optics. 
In  the  following  we  do  not  investigate  the  effects  of 
A(k,  x)  further  because  they  are  similar  to  FTS.  In¬ 
stead,  we  focus  on  the  effect  of  the  phase  term. 


3.  Measuring  Phase  Distortion 

A.  Method 

The  overall  approach  to  determine  phase  distortion  is 
to  use  measurements  of  a  monochromatic  source. 
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Table  1.  Key  Design  Parameters  of  the  SHIMMER-MIDDECK  Instrument 


Attribute 

Design  Value 

Spectral  filter 

Resolving  power 

Spectral  resolution 

Maximum  etendue  (without  optics  transmission) 

2.3  nm  FWHM  centered  at  308.9  nm 

53,500 

58  x  10~3  A 

96  X  10~3  cm2  sr 

For  a  monochromatic  source,  the  modulated  portion 
of  the  interferogram  I  can  be  written  as 

I k0(x)  =  A(k0,  x)(3(K0)(exp{f[K0x  +  <F(k0,  x)]} 

+  exp{-j[K0x  +  <I>(ko,  x)]}).  (2) 

If  the  phase-distortion  term  varies  slowly  with  x,  the 
Fourier  transform  of  I  ,  F (J  ),  has  nonzero  values 
only  in  the  vicinity  of  k0  and  -k0,  corresponding  to 
the  two  exponential  functions  in  Eq.  (2).  By  multi¬ 
plying  F(/  )  with  a  window  function  that  removes  the 
signal  at  and  around  —  k0,  we  can  isolate  the  signal 
around  k0,  and,  after  backward  transformation,  we 
get 

/Ko(x)  =  A(k0,  x)P(K0)exp{i[K0x  +  <b(Ko,  x)]}.  (3) 

The  total  phase  for  k0  can  then  be  obtained  by 

k0x  +  d>(K0,  x)  =  arctan{3[/Ko(x)]/lR[/Ko(x)]},  (4) 

where  2f  and  SR  extract  the  imaginary  and  real  parts 
of  their  respective  arguments.  By  subtracting  the 
first  term  from  the  left  side  of  Eq.  (4),  one  can  deter¬ 
mine  the  phase  distortion  for  the  wave  number  of  the 
monochromatic  source.  To  determine  the  phase  dis¬ 
tortion  as  a  function  of  wave  number,  4>(k,  x),  one  can 
repeat  the  procedure  with  monochromatic  sources  of 
different  wave  numbers.  Most  convenient  is  a  tun¬ 
able  monochromatic  source.  As  an  alternative,  one 
can  use  a  source  with  well-separated  unresolved 
lines.  Isolating  the  different  lines  in  F (/),  as  de¬ 
scribed  above,  then  yields  the  same  result  as  with  a 
true  monochromatic  source. 

At  this  point,  we  also  mention  that  the  shape  of  the 
envelope  function  A(k0,  x)  can  be  obtained  from  Eq. 
(3)  in  that  (3(k0)  is  a  scalar  and  the  exponent  of  the 
exponential  function  is  known  from  the  first  term  in 
Eq.  (4).  The  investigation  of  the  wave-number- 
dependent  envelope  function  is  useful  to  understand 
the  wave-number-dependent  instrumental  line- 
shape  function  of  a  given  interferometer. 

B.  Example 

As  an  example  for  a  measurement  of  the  phase  dis¬ 
tortion  of  a  SHS  spectrometer,  we  use  the 
SHIMMER-MIDDECK  instrument  that  flew  on  the 
Space  Shuttle  Atlantis  (STS- 112)  in  October  20023-5 
and  a  manganese-neon  (MnNe)  hollow-cathode  lamp 
as  the  source.  The  SHIMMER-MIDDECK  instru¬ 
ment  measures  double-sided  interferograms  with  a 
UV-sensitive  CCD.  Some  key  design  specifications 
of  this  instrument  are  listed  in  Table  1.  Figure  2 
shows  the  power  spectrum  (relative  intensity  versus 


fringe  frequency)  of  the  MnNe  lamp  within  the  pass- 
band  of  the  spectrometer  around  309  nm.  We  deter¬ 
mined  the  spectrum  by  transforming  the  modulated 
part  of  the  measured  interferogram  after  correcting  it 
for  detector  offset  (bias  and  dark)  and  a  flat  field. 
From  the  32  interferograms  that  SHIMMER- 
MIDDECK  images  simultaneously  on  the  CCD,3  in¬ 
terferogram  number  16  was  chosen  for  this  example. 
This  interferogram  was  recorded  along  a  CCD  row 
near  the  center  of  the  two-dimensional  array.  The 
power  spectrum  in  Fig.  2  also  includes  the  negative 
spatial  frequencies,  which  result  directly  from  the 
Fourier  transform  of  the  interferogram.  The  spec¬ 
tral  filter  blocks  any  real  signal  at  negative  frequen¬ 
cies  (see  Table  1).  The  spectral  resolution  of  the 
measurement  is  approximately  58  mA. 

The  MnNe  spectrum  features  several  atomic  lines. 
Four  of  these  lines,  identified  as  A,  B,  C,  and  D  in  Fig. 
2,  were  chosen  for  the  determination  of  the  wave¬ 
number-dependent  phase  distortion.  They  were  se¬ 
lected  because  they  are  strong  and  interference  of 
neighboring  lines  is  minimal.  Details  about  the  se¬ 
lected  lines  are  given  in  Table  2. 

Following  the  method  described  above,  we  now  iso¬ 
late  the  individual  lines.  We  chose  to  multiply  the 
complex  spectrum  with  normalized  Gaussian  func¬ 
tions  centered  on  the  maximum  intensity  of  the  lines 
and  a  full  width  at  half-maximum  (FWHM)  of  9.7 
spectral  elements.  The  isolation  functions  are 
shown  in  gray  and  referenced  to  the  right  ordinate  in 
Fig.  2.  Backward  transformation  of  the  isolated 
lines  yields  /  (x)  for  each  line,  and  we  can  now  cal- 


-400  -200  0  200  400 

Spatial  Frequency  (1/detector  array  width) 

Fig.  2.  Black  curve,  power  spectrum  of  a  MnNe  hollow-cathode 
lamp  as  measured  by  the  SHIMMER-MIDDECK  instrument  ver¬ 
sus  spatial  frequency  of  the  recorded  fringes.  Gray  curve,  Gauss¬ 
ian  functions  (FWHM  =  9.7  spectral  elements)  used  to  isolate 
individual  emission  lines  in  order  to  determine  the  frequency- 
dependent  phase  distortion. 
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Table  2.  Emission  Lines  Selected  for  This  Study 


Identifier 

Emitter 

Line  position6  (nm) 

A 

Mn  I 

307.9642 

B 

Ne  II 

308.8170 

C 

Ne  II 

309.2901 

D 

Ne  II 

309.7133 

culate  the  total  phase  [<1>(k0,  x)  +  k0x]  for  each  line  by 
using  Eq.  (4). 

We  define  phase  shift  as  the  phase  at  x  =  0  (near 
the  ZPD).  This  phase  shift  can  be  considered  to  be 
the  value  of  the  phase-distortion  function  at  x  =  0. 
The  phase  shift  for  our  example  spectrum  is  shown  in 
Fig.  3  as  a  function  of  the  wavelength.  The  near- 
linear  wavelength  dependence  is  familiar  from  FTS 
and  is  the  result  of  effects  such  as  the  dispersion 
within  the  interferometer  or  the  sampling  grid,  which 
generally  does  not  include  the  ZPD  location.  For 
increased  clarity  of  the  following  considerations,  we 
subtract  the  phase  shift  from  the  total  phase  of  the 
emission-line  interferograms  so  that  all  measured 
phase  functions  are  zero  at  x  =  0. 

To  determine  the  phase  distortion,  <I>(k0,  x),  from 
Eq.  (4),  we  need  to  determine  k(o-),  where  cr  is  the 
wave  number  of  the  signal.  This  is  equivalent  to  the 
frequency  calibration  of  the  spectrometer.  In  our 
example,  we  fit  a  phase  plane  D[k(o-),  x]  to  the  phase 
functions  of  the  four  emission  lines.  D[k((t),  x]  as¬ 
sumes  a  linear  relationship  between  the  heterodyned 
wave  number  and  k: 

D(ct,  x)  =  k(ct)x  =  C(n  —  cr0)x.  (5) 

The  two  fit  parameters  are  ct0,  the  Littrow  wave  num¬ 
ber,  and  C,  a  constant.  The  result  of  the  fit  is  shown 
in  Fig.  4.  In  this  three-dimensional  plot  the  mea¬ 
sured,  phase-shift-corrected  phase  functions  of  the 
four  emission  lines  are  plotted  as  bold  solid  lines  with 
the  fitted  phase  plane  overplotted  in  gray.  For  the 
Littrow  wave  number  <t0,  the  phase  function  is  zero 
for  each  pixel.  For  pixel  number  0,  near  the  ZPD, 
the  phase  function  is  zero  for  all  wave  numbers  be¬ 
cause  we  phase  shift  corrected  the  measured  phase 
functions.  The  fit  results  in  a  Littrow  wave  number 
o-q  of  32570.4  cm-1  (\  =  307.027  nm)  and  a  C  of 


Fig.  3.  Phase  shift  near  the  ZPD  point  for  the  four  MnNe  emission 
lines.  The  naming  of  the  lines  (A,  B,  C,  D)  follows  from  Fig.  2  and 
Table  2. 


Fig.  4.  Bold,  solid  traces  are  the  phase-shift-corrected  total  phase 
functions  of  the  four  MnNe  emission  lines  listed  in  Table  2.  The 
gray  mesh  is  the  phase  plane  that  is  fitted  to  the  four  phase 
functions.  The  two  fit  parameters  are  the  Littrow  wave  number 
cr0,  where  the  phase  function  is  zero,  and  the  parameter  C,  which 
is  the  linear  slope  increase  of  the  phase  functions  with  respect  to 
the  heterodyned  wave  number  (o-0  —  o);  see  also  Eq.  (4). 


0.0101412  [(rad/cm  ^j/pixel  width],  which  corre¬ 
sponds  to  a  grid  spacing  of  approximately  58  mA  in 
the  spectrum  and  a  maximum  bandwidth  of  approx¬ 
imately  2.97  nm  for  1024  interferogram  samples. 

At  this  point  we  can  determine  the  phase  distortion 
<&(k0,  x),  which,  according  to  Eq.  (4),  is  the  residual  of 
the  measured  phase  functions  and  the  phase  plane 
plus  the  phase  shift  that  we  subtracted  before  the  fit. 

Figure  5  shows  the  phase-shift-corrected  phase  dis¬ 
tortion  for  the  four  emission  lines  in  the  MnNe  spec¬ 
trum.  The  data  show  that  there  is  little  phase 
distortion  around  the  ZPD  location  (x  =  0),  but  a 
small  additional  phase  component  is  apparent  for  all 
four  wave  numbers  toward  higher  optical  path  differ¬ 
ence.  The  data  do  not  show  a  strong  wave-number 
dependence  of  the  phase  distortion,  as  indicated  by 
the  similarity  of  the  four  curves.  The  slight  wave- 


Fig.  5.  Phase-shift-corrected  phase  distortion  for  the  four  emis¬ 
sion  lines  listed  in  Table  2.  All  graphs  are  zero  for  pixel  0  owing 
to  the  phase-shift  correction.  The  phase  distortion  does  not  show 
a  significant  frequency  dependence  (the  curves  basically  fall  on  top 
of  one  another).  The  smoothness  of  the  curves  is  caused  by  the 
narrow  Gaussian  isolation  functions  shown  in  Fig.  2.  The  phase 
distortion  is  small  in  the  center  of  the  detector  array,  and  the 
maximum  phase  distortion  is  less  than  a  fringe  (2ir),  even  for  line 
D,  which  has  the  highest  spatial  frequency  (—467  fringes  per  de¬ 
tector  width)  of  the  four  lines. 
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number  dependence  is  possibly  caused  by  the  inter¬ 
ference  of  neighboring  lines. 

The  smoothness  of  the  phase  distortion  in  Fig.  5  is 
a  direct  consequence  of  the  narrow  Gaussian  isola¬ 
tion  function.  Ideally,  a  monochromatic  source  is 
used  to  determine  the  phase  distortion,  so  that  a 
narrow  isolation  function  is  not  necessary  to  suppress 
neighboring  lines;  thus  the  phase  distortion  can  be 
determined  with  higher  resolution.  However,  as 
long  as  the  phase  distortion  is  sufficiently  slowly 
varying,  a  high-resolution  measurement  of  its  varia¬ 
tion  should  not  be  necessary. 

If  significant  phase  distortion  is  present  for  a  given 
SHS  instrument,  it  is  characteristic  of  the  spectrom¬ 
eter  and  can  generally  be  corrected  for  in  the  post¬ 
processing  of  the  interferogram  with  little  impact  on 
the  signal-to-noise  ratio.  In  the  following,  we  dis¬ 
cuss  several  types  of  phase  error  and  possible  correc¬ 
tion  techniques. 

4.  Phase  Shift  [<t>  =  <|>k(k)] 

The  case  in  which  the  phase  term  in  Eq.  (1)  is  a 
function  of  wave  number  only  [d>  =  ds K  (  k  )  j  is  com¬ 
monly  known  as  phase  shift.  Possible  causes  for 
phase  shift  are  an  uncompensated  beam  splitter  or  a 
sampling  grid  that  does  not  include  the  exact  zero 
path  location  of  the  measured  interferogram.  Phase 
shift  has  been  discussed  widely  in  the  FTS  literature, 
so  we  mention  only  some  key  aspects  of  it  and  refer 
the  reader  to  the  references  for  details. 

Equation  (1)  shows  that  a  phase  shift  is  equivalent 
to  a  multiplication  of  the  recovered  spectrum  by 
exp[i4>K(k)]  for  k  >  0  and  exp[-i4>K(K)]  for  k  <  0.  The 
simplest  correction  method  is  to  calculate  the  absolute 
value  of  the  retrieved  spectrum,  also  called  the  power 
spectrum.  This  operation  ensures  that  the  entire 
spectral  information  is  again  contained  in  the  real  part 
of  the  recovered  spectrum;  however,  it  has  a  nonlinear 
effect  on  the  measurement  noise7  because  both  the  real 
and  the  imaginary  noise  components  contribute  to  the 
power  spectrum. 

Better  correction  techniques  are  based  on  a  deter¬ 
mination  of  4>k(k)  from  the  double-sided  part  of  the 
interferogram  by  use  of  the  fact  that  the  spectrum  by 
definition  has  no  imaginary  part.  Generally,  the  mul¬ 
tiplication  of  the  spectrum8  with  exp[jcj>K(i<)]  or  the 
equivalent  convolution  in  the  interferogram  domain7 
will  yield  the  phase-corrected  spectrum.  Especially 
for  one-sided  interferograms  and  emission-line  spec¬ 
tra,  the  phase-shift  correction  has  to  be  performed 
carefully  to  avoid  the  introduction  of  intensity  and  ab¬ 
solute  line  position  errors.  Further  details  about 
phase-shift  correction  can  be  found  in  the  FTS 
literature.7-10 

5.  Frequency-Independent  Phase  Distortion  [<t>  =  <j>x(x)| 
The  case  in  which  the  phase  term  in  Eq.  (1)  is  a 
function  of  optical  path  difference  only  [T  =  <\>Jx)]  is 
not  typical  for  FTS  but  can  be  found  in  narrow- 
bandpass  SHS  systems.  It  can  be  caused,  for  exam¬ 
ple,  by  grating  figure  errors  or  index  inhomogeneities 
in  the  interferometer.  The  phase  error  <bjx)  can  be 


determined  from  a  single  interferogram  of  one  mono¬ 
chromatic  source  as  described  by  Harlander  et  al2  or 
from  a  multiline  source  as  discussed  previously  in 
this  paper.  Correcting  a  frequency-independent 
phase  distortion2  can  be  achieved  by  a  multiplication 
of  the  interferogram  in  Eq.  (1)  with  exp[— i4>x(x)], 
which  yields 

r(x)^"A(K,x)WK)exp[i(Kx)]+A(K.x,WK) 

X  exp{— j[kx  +  24>3.(x)]}dK.  (6) 

The  first  term  is  identical  to  the  corresponding  term 
in  the  interferogram  without  any  phase  error.  Its 
Fourier  transform  is  the  spectrum  (3(k)  for  k  >  0. 
Here  the  instrumental  line-shape  function  is  the  Fou¬ 
rier  transform  of  the  envelope  function  A(k,  x),  which 
remains  unchanged.  The  second  term  in  Eq.  (6)  ac¬ 
quires  twice  the  phase  error,  and  the  Fourier  trans¬ 
form  of  this  term  results  in  |3(k)  for  k  <  0  convoluted 
with  the  Fourier  transform  of  A(k,  xjexpL— 2id)1.(xj]. 
If  this  corrupted  envelope  function  is  a  slowly  varying 
function  of  x,  its  Fourier  transform  will  be  localized, 
i.e.,  confined  to  low-frequency  components.  In  this 
case  the  second  term  in  Eq.  (6)  will  have  only  a  small 
effect  on  the  k  >  0  part  of  the  Fourier  transform  of 
/'(x),  which  then  contains  the  complete,  corrected 
spectral  information.  The  k  <  0  part  can  be  ignored. 

According  to  the  convolution  theorem  of  Fourier 
transformation,  the  multiplication  with  exp[— i$x(x)] 
in  Eq.  (6)  is  equivalent  to  the  convolution  with  the 
Fourier  transform  of  expf-f^Cx)]  in  the  spectral  do¬ 
main.  For  slowly  varying  phase  distortions  the  con¬ 
volution  kernel  Ffexpf-iclx^x)]}  has  the  advantage  of 
being  localized  in  the  spectral  domain  (see  below). 

6.  Frequency-Dependent  Phase  Distortion  [<t>  =  <t>ix,  k)] 

The  case  of  the  frequency-dependent  phase  distortion 
is  the  most  general  case.  The  previously  discussed 
cases  are  included  in  this  case. 

A  possible  cause  for  frequency-dependent  phase 
distortion  is  image  distortion  by  the  exit  optics  (L2 
and  L3  in  Fig.  1).  If  image  distortion  is  present,  the 
fringe  pattern  at  the  detector  still  appears  in  focus, 
but  the  magnification  of  the  image  changes  across  the 
detector  array.  As  a  consequence,  the  fringe  fre¬ 
quency  from  a  monochromatic  source  changes  across 
the  detector  array.  In  this  case,  the  resulting  phase 
distortion  is  proportional  to  the  fringe  frequency.  In 
our  previously  discussed  example,  image  distortion 
would  cause  the  phase  distortion  of  line  D  (—467 
fringes  per  detector  width)  to  be  approximately  2.8 
times  the  phase  distortion  of  line  A  (—164  fringes  per 
detector  width).  Figure  5  does  not  show  this  behav¬ 
ior,  which  indicates  that,  for  this  instrument,  image 
distortion  is  not  the  dominant  contribution  to  the 
phase  distortion. 

The  correction  of  the  frequency-dependent  phase 
distortion  has  to  be  performed  in  the  spectral  domain 
because  the  correction  is  different  for  each  wave  num¬ 
ber  and  the  interferogram  samples  carry  information 
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of  all  wave  numbers  at  the  same  time.  In  the  case  of 
the  frequency-independent  phase  distortion  we 
found,  at  the  end  of  Section  5,  that  the  correction  can 
be  performed  in  the  spectral  domain  via  a  convolution 
with  a  correction  function.  Here  we  will  apply  a 
similar  technique. 

To  correct  for  the  frequency-dependent  phase  dis¬ 
tortion,  we  first  calculate  the  Fourier  transform  of 
exp[-i4>(x,  k)]  for  each  k,  which  we  call  Gk(k'): 

Gk(k')  =  F(exp{-i[4>(x,  k)]}).  (7) 

Then  we  correct  the  spectrum  [3  r/(  k  ' )  that  results 
from  the  Fourier  transformation  of  the  uncorrected 
interferogram  by  convolving  it  with  Gk(k'): 

Pc(k)  =  J  Mk")Gk(k  -  K")dK".  (8) 

In  the  strict  mathematical  sense,  this  is  only  a  con¬ 
volution  if  the  function  Gk(k')  is  independent  of  k, 
which  is  the  case  of  frequency-independent  phase  dis¬ 
tortion,  and  it  is  equivalent  to  the  correction  tech¬ 
nique  previously  discussed  for  that  case.  Here  the 
convolution  function  depends  on  k.  However,  if  the 
phase  distortion  is  slowly  varying  with  x  and  k,  Gk(k') 
will  also  vary  slowly  with  k,  and  all  Gk(k')  functions 
are  nonzero  only  around  k'  =  0.  This  means  that 
locally  (around  each  k  in  the  spectrum)  the  conditions 
are  similar  to  a  real  convolution,  and  good  correction 
results  can  be  achieved  with  this  technique.  In  Sec¬ 
tion  7  we  will  illustrate  this  technique  by  using  the 
highly  structured,  near-UV  solar  spectrum  as  an  ex¬ 
ample. 

7.  Simulated  Performance  of  Frequency-Dependent 
Phase-Distortion  Correction  with  a  Near-UV  Solar 
Spectrum 

To  assess  the  performance  of  the  frequency- 
dependent  phase-distortion  correction  discussed 
above,  we  take  a  high-resolution  near-UV  solar  spec¬ 
trum,11  multiply  it  with  a  filter  transmittance  cen¬ 
tered  at  approximately  309  nm,  and  bin  it  into  513 
frequency  bins  with  a  width  of  60  mA.  These  pa¬ 
rameters  have  been  chosen  because  they  are  similar 
to  the  parameters  of  the  SHIMMER  instruments.3-12 
By  use  of  Eq.  (2)  as  a  model,  synthetic  interferograms 
were  generated  for  this  spectrum  for  the  various 
phase-distortion  functions  described  below.  Figure 
6,  panel  A  (6A),  shows  the  spectrum  without  any 
phase  distortion. 

The  first  simulation  uses  a  set  of  phase-distortion 
functions  that  are  a  superposition  of  a  frequency- 
independent  contribution  and  a  function  that  is  pro¬ 
portional  to  the  frequency  and  to  the  fourth  power  of 
the  pixel  number  to  simulate  image  distortion.  Fig¬ 
ure  6B  shows  three  phase-distortion  functions  for  the 
frequencies  marked  o-j,  an,  and  o-m  in  Fig.  6A.  The 
leftmost,  black  marking  in  Fig.  6A  corresponds  to  the 
Littrow  frequency.  Figure  6C  shows  the  real  part  of 
the  uncorrected  spectrum  that  is  retrieved  by  the 
Fourier  transform  of  the  phase-distorted  interfero- 
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Fig.  6.  First  simulation  case:  high  spectral  resolution  and  non¬ 
periodic  phase  distortion.  A:  high-resolution  solar  spectrum 
multiplied  with  a  typical  interference  filter  transmittance.  The 
leftmost,  black  mark  indicates  the  Littrow  frequency,  and  gray 
marks  indicate  the  frequencies  for  which  the  phase-distortion  func¬ 
tions  are  plotted  in  panel  B.  B:  phase-distortion  functions  for 
the  frequencies  marked  in  panel  A.  C:  black  curve,  real  part  of 
the  uncorrected,  phase-distorted  spectrum.  Gray,  residuals  be¬ 
tween  the  uncorrected  spectrum  and  the  initial  spectrum.  D: 
black  curve,  phase-distortion-corrected  spectrum.  Gray,  residu¬ 
als  between  the  phase-distortion-corrected  spectrum  and  the  ini¬ 
tial  spectrum  multiplied  by  10. 


gram.  It  also  shows  the  residuals  with  respect  to  the 
initial  spectrum  (Fig.  6A)  in  gray.  Figure  6D  shows 
the  phase-corrected  spectrum  and  the  residuals  mul¬ 
tiplied  by  10  in  gray.  The  correction  has  improved 
the  residuals  by  more  than  an  order  of  magnitude, 
but  a  systematic  residual  at  the  1%  level  is  still  ap¬ 
parent.  The  correction  is  not  complete  because  it 
relies  on  the  assumption  that  the  convolution  func¬ 
tions  (the  Fourier  transform  of  the  phase-distortion 
functions)  are  slowly  varying  with  frequency  and 
quickly  falling  off  to  zero  away  from  the  center.  The 
phase-distortion  functions  and  corresponding  convo¬ 
lution  functions  that  we  used  for  the  first  simulation 
violate  mainly  the  second  assumption.  The  convo¬ 
lution  functions  in  this  simulation  do  not  show  rapid 
convergence  to  zero  away  from  the  center  because  the 
phase  distortion  for  the  lowest  frequency  is  signifi¬ 
cantly  different  from  that  at  the  highest  frequency 
and  the  discrete  Fourier  transform  of  this  disconti¬ 
nuity  results  in  significant  ringing. 

One  can  avoid  this  behavior  by  choosing  phase- 
distortion  functions  that  have  the  same  value  for  the 
highest  and  lowest  frequencies,  as  seen  in  Fig.  7B. 
As  a  result,  the  residuals  (Fig.  7D)  of  the  corrected 
spectrum  improve  significantly.  Of  course,  this  is 
not  a  real  option  for  a  given  spectrometer,  in  which 
the  phase-distortion  functions  cannot  be  changed 
easily.  Another,  more  practical  way  to  improve  the 
performance  of  the  correction  is  to  reduce  the  resolu¬ 
tion  of  the  spectrum,  for  example  by  apodizing  the 
interferogram.  The  lower-resolution  spectrum  will 
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Fig.  7.  Second  simulation  case:  high  spectral  resolution  and  pe¬ 
riodic  phase  distortion.  A:  high-resolution  solar  spectrum  mul¬ 
tiplied  with  a  typical  interference  filter  transmittance.  The 
leftmost,  black  mark  indicates  the  Littrow  frequency,  and  the  gray 
marks  indicate  the  frequencies  for  which  the  phase-distortion  func¬ 
tions  are  plotted  in  panel  B.  B:  phase-distortion  functions  for 
the  frequencies  marked  in  panel  A.  C:  black  curve,  real  part 
of  the  uncorrected,  phase-distorted  spectrum.  Gray,  residuals 
between  the  uncorrected  spectrum  and  the  initial  spectrum.  D: 
black  curve,  phase-distortion-corrected  spectrum.  Gray,  resid¬ 
uals  between  the  phase-distortion-corrected  spectrum  and  the 
initial  spectrum  multiplied  by  100.  The  residuals  are  more 
than  an  order  of  magnitude  smaller  than  in  the  first  simulation 
case  (Fig.  6). 


effectively  minimize  the  contributions  of  the  high- 
frequency  ringing  of  the  convolution  function  in  Eq. 
(8).  However,  the  price  of  lower  spectral  resolution 
has  to  be  paid.  Figure  8A  shows  the  solar  spectrum 
convoluted  with  the  line-shape  function  that  corre¬ 
sponds  to  an  interferogram  apodization  with  a  Han¬ 
ning  function.  For  this  case,  the  residuals  of  the 
corrected  spectrum  are  well  below  the  1%  level  (see 
Fig.  8D). 

Many  spectra  are  calculated  with  apodized  inter- 
ferograms  from  the  start  to  achieve  a  well-localized 
instrumental  line-shape  function  rather  than  a  sine 
line-shape  function  that  also  includes  significant 
ringing.  In  case  the  loss  in  resolution  by  apodizing 
the  interferogram  is  not  acceptable,  the  phase- 
distortion  term  can  be  included  in  a  frequency- 
dependent  instrumental  line-shape  function  for 
further  data  analysis. 

8.  Conclusion 

Interferograms  measured  by  SHS  instruments  may 
show  phase  distortion  depending  on  the  instrument 
design  and  implementation.  If  significant  phase 
distortion  is  present  for  a  given  instrument,  it  can 
generally  be  corrected  without  significant  loss  of  the 
signal-to-noise  ratio.  We  presented  a  method  to  de¬ 
termine  the  frequency-dependent  phase  distortion 
and  applied  it  to  a  measurement  of  a  UV  SHS  instru¬ 
ment  by  using  a  multiline  source.  For  this  case,  we 
found  the  phase  distortion  to  be  small  (less  than  a 
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Fig.  8.  Third  simulation  case:  lower  spectral  resolution  and 
nonperiodic  phase  distortion.  A:  high-resolution  solar  spectrum 
convolved  with  the  kernel  (i,  i,  1),  which  is  equivalent  to  an  inter¬ 
ferogram  apodization  with  a  Hanning  function,  multiplied  with  a 
typical  interference  filter  transmittance.  The  leftmost,  black 
mark  indicates  the  Littrow  frequency,  and  the  gray  marks  indicate 
the  frequencies  for  which  the  phase-distortion  functions  are  plot¬ 
ted  in  panel  B.  B:  phase-distortion  functions  for  the  frequencies 
marked  in  panel  A.  C:  black  curve,  real  part  of  the  uncorrected, 
phase-distorted  spectrum.  Gray,  residuals  between  the  uncor¬ 
rected  spectrum  and  the  initial  spectrum  multiplied  by  10.  D: 
black  curve,  phase-distortion-corrected  spectrum.  Gray,  residu¬ 
als  between  the  phase-distortion-corrected  spectrum  and  the  ini¬ 
tial  spectrum  multiplied  by  100.  The  residuals  are  more  than  an 
order  of  magnitude  smaller  than  in  the  first  simulation  case  (Fig. 
6)  and  well  below  the  1%  level. 


fringe  even  for  the  highest  spatial  frequencies)  and 
virtually  frequency  independent.  We  subsequently 
simulated  three  phase-correction  cases  by  using  the 
UV  solar  spectrum  around  309  nm  with  different 
spectral  resolutions  and  different  phase-distortion 
functions.  We  verified  that  the  correction  method 
works  well  if  the  phase  distortion  varies  slowly  with 
frequency  and  if  the  Fourier  transform  of  the  phase- 
distortion  functions  falls  off  quickly  for  higher  fre¬ 
quencies.  Since  the  second  assumption  is  not 
necessarily  true,  e.g.,  if  the  phase-distortion  func¬ 
tions  are  not  periodic,  we  showed  that  decreasing  the 
resolution  of  the  spectrum  (apodizing  the  interfero- 
gram)  improves  the  quality  of  the  correction. 

Small  phase  distortions  like  the  ones  in  our  exam¬ 
ple  should  be  the  goal  of  every  SHS  spectrometer 
design  but  are  not  always  easy  to  achieve.4  In  these 
cases  the  proper  correction  gains  in  importance. 
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